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1 Introduction 

Let Xi , . . . , Xm ,Yi, . . . ,Yn be independent non- negative random variables such that each Xi 
has the probability distribution Pi and each Yj has the probability distribution P2. Given 
realized values X = {Xi}^^ and Y = {Yj}'j^^ we define the random bipartite graph Hx,y 
with the bipartition V = {vi, . . . ,Vn}, W = {wi, . . . jWm}, where edges {wi^vj} are inserted 
with probabilities pij = mm{l, XiYj{n'm)~^^'^} independently for each {i,j} G [m] x [n]. The 
inhomogeneous random intersection graph G{Pi, P2, n, m) defines the adjacency relation on the 
vertex set V: vertices v' ,v" G V are declared adjacent (denoted v' ~ v") whenever v' and v" 
have a common neighbour in Hx,y- 

The degree distribution of the typical vertex of the random graph G{Pi, P2,n, m) has been first 
considered by Shang [20j. The proof of the main result of [20] contains gaps and the result is 
incorrect in the regime where m/n — )• /3 G (0, -|-cxd) as m, n — )• +00. We remark that this regime is 
of particular importance, because it leads to inhomogeneous graphs with the clustering property: 
the clustering coefficient P{vi ~ V2\vi ~ vs,V2 ~ U3) is bounded away from zero provided that 
EXf < 00 and EY^ < 00, see [H]. The aim of the present paper is to show the asymptotic 
degree distribution in the case where m/n — t- /3 for some f3 G (0, +00). 

We consider a sequence of graphs {Gn = G{Pi, P2,n,m)}, where m = rrin +00 as n — ^ 00, 
and where Pi,P2 do not depend on n. We denote = EX{, bi = EiY^. By d{vj) = dG„{vj) 
we denote the degree of a vertex Vj in G„ (the number of vertices adjacent to vj in Gn)- We 
remark that for every n the random variables dG„{'^i)j ■ ■ ■ ■,dQ^{vn) are identically distributed. 
In Theorem [1] below we show the asymptotic distribution of d{vi). 

Theorem 1. Let m,n — t- cxd. 



(i) Assume that m/n — t- 0. Suppose that that EXi < oo. Then P(d(t)i) = 0) = 1 — o(l). 

(ii) Assume that m/n ^ fi for some j3 G (0, +oo). Suppose that EX^ < oo and EYi < oo. Then 
d{vi) converges in distribution to the random variable 



Ai 

(1) 
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where ti,T2, . . . are independent and identically distributed random variables independent of the 
random variable Ai. They are distributed as follows. For r = 0, 1, 2, . . . , we have 

^i^i = r) = '^P{M = r + l) and P(Ai = r) = E e"^' ^, i = l,2. (2) 
EA2 r! 

Here Ai = Yiail3^/^ and A2 = Xibip-^/^. 

(Hi) Assume that m/n — )• +00. Suppose that EXf < 00 and EYi < 00. Then d{vi) converges 
in distribution to a random variable A3 having the probability distribution 

P(A3 = r) = Ee-^3^, r = 0,l,.... (3) 
r! 

Here A3 = Yia2bi. 

Remark 1. The probability distributions of Aj, i = 1,2,3, are Poisson mixtures. One way 
to sample from the distribution P\. is to generate random variable Xt and then, given Aj, to 
generate Poisson random variable with the parameter Aj. The realized value of the Poisson 
random variable has the distribution P{^.. 

Remark 2. The asymptotic distributions ([T]) and (31) admit heavy tails. In the case (ii) we 
obtain a power law asymptotic degree distribution (l| provided that at least one of Pi and P2 
has a power law and Pi(0), ^2(0) < 1. In the case (iii) we obtain a power law asymptotic degree 
distribution ^ provided that P2 has a power law. 

Remark 3. Since the second moment 02 does not show up in ([T]), ([2]) we expect that in the case 
(ii) the second moment condition EXf < 00 is redundant and could perhaps be replaced by the 
weaker first moment condition 'EXi < 00. 

Random intersection graphs have attracted considerable attention in the recent literature, see, 
e.g., [5, [2], [3], [9], [TO], [12], [E]. Starting with the paper by Karohski et al [TO], see also [2T] . 
where the case of degenerate distributions Pi = Pi„, P2 = P2n depending on n was considered 
(i.e., Pin(cn) = P2n(cn) = 1, for some Cn > 0), several more complex random intersection graph 
models were later introduced by Godehardt and Jaworski [13], Spirakis et al. [T7j, Shang [20]. 
The asymptotic degree distribution for various random intersection graph models was shown in 

m, m, m, m, m, m, m, m, m- 

2 Proofs 

Before the proof we introduce some notation and give two auxiliary lemmas. 

The event that the edge {wi, Vj} is present in H = Hx,y is denoted Wi — )• Vj. We denote 

m 

hj=\w,^Vj}, Ii=Iil, = hv L = ^^Uili. 

2<j<n i=l 
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We remark, that Ui counts all neighbours of Wi in H belonging to the set V \ {vi}. Denote 



ak = m-^ Xl h = n-^ ^ y/, 

l<.i<.m 2<j<n 

Xij = SxY = VI -i=^^i3 min{l, Aij}, (4) 

and introduce the event A\ = {Xn < 1, 1 < i < m}. By Pi and Ei we denote the condi- 
tional probability and conditional expectation given Yi. By P and E we denote the conditional 
probability and conditional expectation given X, Y. By dxyiC^ C) we denote the total variantion 
distance between the probability distributions of random variables Q and ^. In the case where ^ 
and X^Y are defined on the same probability space we denote by dTviCiCj the total variation 
distance between the conditional distributions of C and ^ given X, Y . 

In the proof below we use the following simple fact about the convergence of a sequence of 
random variables 

Xn = op{l), Esupx„<oo =^ Ex„ = o(l). (5) 

We remark that the condition E sup„ x„ < oo (which assumes implicitly that all Xn are defined 
on the same probability space) can be replaced by the more resctrictive condition that there 
exists a constant c > such that < c, for all n. In the latter case we do not need all x„ to 
be defined on the same probability space. In particular, given a sequence of bivariate random 
vectors {{rjniOn)} such that, for every n and m = run random variables r/„, 0n and {Xi}'^^, 
{Yj}^^^ are defined on the same probability space, we have 

dTv{r]n,On) = Op{l) dTviVn^On) <EdTv{'nn,On) = 0{l). 

Lemma 1. Assume that EX^ < oo and EYi < oo. We have as n,m ^ +oo 

P{d{vi) / L) = 0(1), (6) 
P(A) = 1 - 0(1), (7) 
SxY = op{l), ESxY = o{l). (8) 

Proof of Lemma^ Proof of ([6]). We observe that the event A = {d{vi) ^ L\ occurs in 
the case where for some 2 < j < n and some distinct ii,i2 G [m] the event Aij^^i2j = 
{wi-^ — )• vi, Wi^ — )• Vj, Wi^ — )■ fi, Wi^ — )• Vj} occurs. From the union bound and the inequal- 
ity P(Ai,i2,i) < m-^n-^XlXlY^Yf we obtain 

P(^) = P ( U U I < n-%Y,^Qx. (9) 

\{n,i2}cM2<j<n / 

Here Qx = m^'^ Yl{ii i2}c[m] -^h-^lz' note that Qx and Y^ are stochastically bounded and 
n~^b2 = op(l) as n — )• +oo. Therefore, P(^) = op(l). Now implies wm. 



Proof of ([7j). Let Ai denote the complement event to Ai- By the union bound and Markov's 
inequality 

Pi(:4i) < ^ 1) ^ {nmy^Y^ ^ = ^"'«2n'- 

iG[ml is [ml 



3 



Hence we obtain Pi(A) = op(l). Now ^ implies P(A) = EPi(^i) = o(l). 

Proof of Q. Since the first bound of ([s]) follows from the second one, we only prove the latter. 
Denote Xi = max{Xi, 1} and Yi = max{Yi, 1}. We observe that EXf < oo, EYi < oo implies 

lim BXhri> = 0, lim BYiIr^ = 0. 

Hence one can find a strictly increasing function ip : [l,+oo) — )• [0,+oo) with limt^^oo (p{t) = 
+00 such that 

EXiV(Xi) < oo, Eyi^(li) < oo. (10) 
In addition, we can choose ip satisfying 

ip{t)<t and ip{st) < Lp{t)ip{s), Vs,t>l. (11) 

For this purpose we take a sufficiently slowly growing concave function ^ : [0, +oo) — t- [0, +oo) 



with V(0) = and define (p{t) = e'^^'"^*^^ We note that the second inequality of (11) follows 



from the concavity property of tp. We remark, that (11) implies 



t/{st) = s''^ <l/<f{s) <<f{t)/Lp{st), s,t>l. (12) 



Let SxY be defined as in (|4j) above, but with Ajj replaced by Ay = XiYj / y/rrm. We note, that 
SxY < SxY- Furthermore, from the inequalities 

„i„0,A,,a<miJl,4p>}<^<?^(|M5^ (13, 

we obtain Sxv < Sxy < S^y/(^(\/nm), where 

S*xY = {mn)-A XfifixA I YMYj. 

\l<i<m I \2<j<n 



We remark that (11) and (12) and imply the third and the first inequality of (13), respectively. 



Finally, the bound ESxy = o(l) follows from the inequality Sxy < / fiV^^^) ^^^^ the fact 



that ES'^y remains bounded as n, m — )• +oo, see (10). □ 



In the proof of Theorem [T] we use the following inequality refered to as LeCam's lemma, see e.g. 



Lemma 2. Let S = Ii + I2 + • • • + In be the sum of independent random indicators with 
probabilities P(Ii = 1) = pi. Let A be Poisson random variable with mean pi + ■ ■ ■ + Pn- The 
total variation distance between the distributions Pg of Pa of S and A 

sup |P(5 G A) - P(A G ^)| = J V |P(5 = k)-V{k = k)\< (14) 

AC{0,1,2...} . 

Proof of Theorem^ The case (i). We have F{d{vi) > e) < P(-^ > £), for each e > 0. We prove 
that P(L > e) = o(l), for any e > 0. In view of the identity P(L > e) = EPi(L > e) and ^ 
it suffices to show that Pi{L > e) = op(l). For this purpose we write, by the union bound and 
Markov's inequality, 

Pi{L>e)< J2 Pi(Ifc = l)<Ei Y >^ki = VWnYi-EXi = op{l). 

l<fc<m l<fc<m 
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In cases (ii) and (iii) we apply ([6]). In view of ([g]) the random variables d{vi) and L have the 
same asymptotic distribution (if any). Hence, it suffices to show the convergence in distribution 
of L. 

The case (ii). Here we prove that L converges in distribution to ([T]). We first approximate L 
by the random variable L3 = Y^kLi "HkCsk- Then we show that L3 converges in distribution to 
([1]). Here r/i, . . . , r/m, ^31, ... , ^3^ are conditionally independent (given X, Y) Poisson random 
variables with Er/jt = A^i and £^3^ = Xfc(n/m)^/^6i. We assume, in addition, that given X,Y, 
the sequences {Ifcj^x ^^'^ {Csfcl^Li conditionally independent. 

Given X, Y, we generate independent Poisson random variables ^n, . . . , ^im, An, . . . , Ai^,, with 
the conditional mean values 

E6a: = ^ Pkj, EAifc = ^ {>^kj-Pkj), l<k<m. 

2<j<n "iKjKn 

We assume that, given X, Y , these Poisson random variables are conditionally independent of the 
sequence We suppose, in addition, that {%}^]^ is conditionally independent (given 

X,Yi) of the set of edges of H that are not incident to vi. We define ^2fc = £,ik + Ai^ and 
observe that ^2fc has conditional (given X, Y) Poisson distribution with the conditional mean 
value E^2fc = ^2<j<n^kj- Introduce the random variables 

l<k<.m l<k<m l<k<m 

In order to show that L and L3 have the same asymptotic probability distribution (if any) we 
prove that 

drviL, Lo) = 0(1), drviLo, L{) = o(l), (15) 
EIL1-L2I =0(1), L2-L3 = op(l). (16) 

Here L2 and L3 are marginals of the random vector [L2,L^) constructed below which has the 
property that L2 has the same distribution as L2 and L3 has the same distribution as L3. 



Let us prove the first bound of (15). In view of (pj) it suffices to show that dTvi^Q, L) — op(l) 



In order to prove the latter bound we apply the inequality 

dTviLo,L)I_A, <n-'Y^ d2 (17) 



shown below. We remark that (17) implies 

dTv{Lo,L) < dTv{Lo,L)lA,+I^^ < n-%^d2+I^^ = op(l). 

Here n~^Yid2 = op(l), because Y^d2 is stochastically bounded. Furthermore, the bound 
op(l) follows from Q. 



It remains to prove (17). We denote L'^ = Yli=i^i'^i + SilLfc+i and write, by the triangle 
inequality, 

m 

dTviLo,L) <^dTv{L'k-i^L'^)- 

k=l 

Then we estimate dTv{L'k_i, L'j^) < drviVkj^k) < {nm)~^YiX^. Here the first inequality follows 
from the properties of the total variation distance. The second inequality follows from Lemma 
[2] and the fact that on the event Ai we have pki = Xki- 
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Let us prove the second bound of ( 15). In view of ([5j) it suffices to show that dTvi^o, Li) = op(l) 
We denote = Yli=i 'Hi'^i + Yl'iLk+i ViCu and write, by the triangle inequahty, 



dTviLo,Li) < y^^dTviLl_i,Ll). 



(18) 



k=l 



Here 



dTv{Ll^i,Ll) < dTv{VkUk,m^ik) < i'iVk 7^ 0)dTv{uk,(ik)- (19) 

Now, invoking the inequahties P^rjk / 0) = 1 — e~^''^ < Xki and dTv{uk,^ik) < Yl^=2Pkj^ 
obtain from ([is]), ([l9]) and ([s]) that 

m n 

k=l j=2 



Let us prove the first bound of (16). We observe that 



\L2 - Li\ = L2 - Li = ^ rjkAik 



l<k<m 



and 



E J2 "^kAik = ^^kj - l)^{Afe,>l} ^ YlSxY. 

l<k<m l<k<m 2<j<n 

We obtain B\L2 - Li\ < BYiESxY = o(l), see 

Let us prove the second bound of ( |16| ). Given X,Y, generate independent Poisson random 
variables ^3;^, . . . , ^3^, A21, . . . , A2m, A31, . . . , A-^m which are conditionally independent of the 
sequence {r]k}k^=i and have the conditional mean values 

E^sfc = Xk{n/my/\ EAafe = Xk{n/m)^/H2, EAgfc = Xk{n/m)^/^63. 

Here b = min{6i,6i}, 82 = bi — b, 5^ = bi — b. We note that 82,53 > and observe that the 
random vector 

(L'2,4), 4= Y r]ki^3k + ^2k), 4= Y %(ek. + A3fe) 

l<k<m l<fe<m 

has the marginal distributions of (L2,L3). In addition, we have 

E|L'2 - 41 = {62 + 53)^102 = |6i - ^ilnas = op(l). (20) 

In the last step we used the fact that Yid2 = Op(l) and 61 — &i = op(l), by the law of large 
numbers. Finally, we show that (20) implies the bound ~ -^3! = op(l). Denoting, for short, 
H = 14 — 41 and h = El^fj^^^y we write, for e £ (0, 1), 



P{H >e)=Eh = E/i(I^E/f>e2} + hEH<e^})- 



(21) 



Using the simple inequality h < 1 and the inequality, h < e ^EH, which follows from Markov's 
inequality, we obtain 

EM{Eff>e2} < EI^E/f>e2} = P(Ei^ > e') = 0(1), 
^h\EH<e^ < E(e-iEi/)I^EH<.2} < e. 
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Invoking these inequalities in (21) we obtain P{H > e) < e + o(l). Hence H = op(l). 
Now we prove that L3 converges in distribution to ([l]). Introduce the random variable L = 
^i<k<m VkCk, where, given X, Y, the random variables ^1, . . . , are conditionally independent 
of {rik}^^i and have the conditional mean values E^fc = Xkf3~^^'^bi. Proceeding as in the proof of 
the second bound of (16) above, we construct a random vector (L3, L') with the same marginals 
as (L3, L) and such that 



L'\ < |1 - {m/n)^/'^p-^/^\Yibia2 = op(l). 



(22) 



In the last step we used the fact that lia2 = Op(l) and m/n — )• /3. Now, (22) implies L3 — L' 



op(l). We conclude that L3 and L have the same asymptotic distribution (if any). 
Next we prove that L converges in distribution to |lj). For this purpose we show that Ee**^ — ?■ 
Ee**''*, for each t G (—00, +00). Denote A(t) = — e**''*. We shall show below that, for any 
real t and any realized value Yi there exists a positive constant c = c{t,Yi) such that for every 
< e < 0.5 we have 

limsup |E(A(t)|yi)| < ce. (23) 



Clearly, (23) implies E(A(t)|yi) = o(l). This fact together with the simple inequality |A(f)| < 2 
yields EA(t) = o(l), by Lebesgue's dominated convergence theorem. Observing that EA(t) = 
jjgiiL _ ^^itd, conclude that Ee**-^ — )- Ee**''* . 



We fix < e < 0.5 and prove (23). Before the proof we introduce some notation. Denote 

m 



Pr 



r>0 l<k<m 

S = {fAt)-l)\-{Ir{t)-l)\u /(t) = E,e''''-, /(() 



k=l 



Ee 



itL 



Here E denotes the conditional expectation given X,Y and ■ ■ ■ ,irn- Introduce the event 
T) = {|di — Oil < eminjl, ai}} and let V denote the complement event. Furthermore, select the 
number T > 1/e such that P(ri > T) < e. By ci, C2, . . . we denote positive numbers which do 
not depend on n, m. 

We observe that, given Yi, the conditional distribution of d^, is the compound Poisson distribu- 
tion with the characteristic function f{t) = e'*'^^^^^^^~^\ Similarly, given X, Y and ^1, . . . , the 
conditional distribution of L is the compound Poisson distribution with the characteristic func- 
tion f{t) = e'^(-^^W~^). In the proof of (23) we exploit the convergence A — >• Ai and frit) — )• frit)- 
In what follows we assume that m, n are so large that /3 < 2m/n < 4/3. 



Let us prove (23). We write 

EiA(t) = /i + /2. 



h = EiA(t)fe, h = EiA(t)%. 



Here I/2I < 2Pi(D) = 2P(D) = o(l), by the law of large numbers. Next we estimate Ii. 
Combining the identity EiA(t) = Eiif{t){e^ - 1) with the inequalities |/(t)| < 1 and |e* - 1| < 
|s|el''l, we obtain 

< Ei\6\e^^kv < ciEil^lIi,. (24) 
Here we estimated e'^' < e®^^ =: ci using the inequalities 



|5|<2A + 2Ai, A = yi(m/n)^/2^i < 3Ai. 
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Wc remark that the last inequality holds for m/n < 2/3 provided that the event T> occurs. 
Finally, we show that Ei|(5|Ix) < (c2 + A1C3 + AiC4)e + o(l). We first write 

S = (frit) - 1)(A - Ai) + (frit) - fr{t))Xl, 

and estimate \6\ < 2|A — Ai| + Ai|/r(t) — Prom the inequalities |ai — ai| < e and m/n < 2/3 
we obtain 

|A - All < Yi\ai - a|(m/n)V2 + Yiai\{m/n)^/^ - < 2Yip^/^e + o(l). 
Hence Ei|A — Ai|Id < C2£ + o(l), where C2 = 2Yi(3^/'^. We secondly show that 

Bl\fr{t) - frimv < (C3 + 04)6 + 

To this aim we split 

frit) - frit) = ^'"iPr " Pr) = Rl - R2 + R3, 
r>0 

and estimate separately the terms 

R, = Ye'i'pr, R2 = Y.e'''Vr, Rs = Yl ^'"^Pr-Vr). 

r>T r>T 0<r<T 

Here we denote Pr = P(ri = r). The upper bound for R2 follows by the choice of T 

\R2\ < =P(ti >r) <£. 

r>T 

Next, combining the identity pr = (fiirn)"^ X^i<fc<m "'^fc^{4=r} with the inequalities 

\Ri\<{aim)-'Yl E Xk\^,=r} = i^irn)-' E ^fe%fc>T} (25) 

r>Tl<k<m l<k<m 

and Ui^Ix) < ^a^^, we estimate 

^i\Ri\lv < 2a^^Ei(XiI|^-^>^|) < 2a^^r-iEi(Xi^i) = 2a^^a2bi/3~^/^£- 
Hence Eiji^illx" < ca^, where C4 = 2a^^a2bi/3^^/'^ . 

Now we estimate R3. We denote p'^ = {ai/ai)pr and observe that the inequality |ai — ai| < eai 
implies |aia]^^ — 1| < e and 

I E ^'''{Pr-P'r)\<e Pr<e. 
0<r<T 0<r<T 

In the last inequality we use the fact that the probabilities {pr}r>o sum up to 1. It follows now 
that 

\R3\lv<£+ Y \Pr-Pr\- 
0<r<T 

Furthermore, observing that Eip^ = a~^'EiXkI^^^^^y = pj., for 1 < k < m, we obtain 
^i\Pr-Prf = m-^'^i\a-^XiI^^^^^^ -pr\'^ < m-^a^'^EXf. 
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Hence, Ei|p^ — Pt\ = 0{m ^/'^). We conclude that 

Ei|i?3|I© <e + 0(|T|m-i/2) ^ ^ ^ ^^-^y 

The case (in). We start with introducing some notation. Denote m/n = /3n. Given e G (0, 1) 
introduce random variables 
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l<k<m 



Given X,Y, let Ii, . . . , Im be conditionally independent Bernoulli random variables with success 
probabilities 

p(ifc = 1) = 1 - p(ifc = 0) = 7fc. 

We assume that, given X,Y, the sequences {Ifej^i, {Ifcl^Li {CsfcjfcLi conditionally 
independent. Introduce random variables 



l<fc<m 



l<A,<m 



l<k<m 



Furthermore, we define the random variable L7 as follows. We first generate X, Y. Then, given 
X,Y, we generate a Poisson random variable with the conditional mean value 7. The realized 
value of the Poisson random variable is denoted L7. Thus, we have F^Lj = r) = 'Eie~^^'^ /r\, for 
r = 0,l,.... 

We note that L and -L3 have the same asymptotic distribution (if any), by (15), (16). Now we 
prove that L3 converges in distribution to A3. For this purpose we show that for any e £ (0, 1) 



dTv{L3, U) = 0(1), E(L4 - L5) = 0(1), 
dTviL^^Lo) < a2b{e, drviLQ^L-j) = o(l) 
Ee**-^^ - Ee**^^' = oil). 



(26) 
(27) 
(28) 



Let us prove (26), (27), (28). The first bound of (26) is obtained in the same way as the first 
bound of (15). To show the second bound of (26) we invoke the inequality 

E(L4-L5)= (1 - I^EIfciEesfc < n^im-i (1-4)^' 



l<fc<m 



l<fc<m 



and obtain 



E(L4 - L5) = EE(L4 - U) < 6?EX2l , 1/2, = o{l). 



We note that the right hand side tends to zero since /3„ — )■ +00. 



Let us prove the first inequality of (27). Proceeding as in (18), (19) and using the identity 
Ik = Ikl'k we write 



dTviL^,LG)< Y IfcP(Ifc/0)(iry(C3fc,lEfc)- 



l<k<m 



Next, we estimate I^fiTy (Csfc; Ife) < 7fc, by LeCam's inequality (14), and invoke the inequality 
P(Ifc / 0) < Afci. We obtain 

dTv{L5,LQ)< Y 4^A:i7fc < e Y 4-^^17^ < £^1^102- 

l<fc<m l<k<m 
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Here we estimated 7I < ejk- Now the inequalities clxviL^, Lq) < 'Edrvi^^, Lq) ^ a2^ie imply 
the first relation of (27). 

Let us prove the second relation of (27). In view of (p| it suffices to show that dTviLe, L^) = 



op(l). For this purpose we write 



where = op(l), see (j7|), and estimate using LeCam's inequality (14) 

1^1 dry (Le, Lj) < 1^, Yl ^^O^kh = 1)4 < 6? if ^-2^4 = op(l) 



l<k<m 



Here we used the fact that EX^ < 00 implies m ^d4 = op(l). 
Finally, we show (28). We write Ee**^^ = e'''*^^'*"^) and observe that 

Yibia2 - 7 = op{l). 



(29) 



Furthermore, since for any real t the function z — )• e^^*^'* is bounded and uniformly continuous 
for 2; > 0, we conclude that (29) implies the convergence 



It remains to prove (29). We write 116102 — 7 = Yibi{a2 — 02) + li6i«2 — 7 and note that 
0^2 — 02 = Op (1) 5 by the law of large numbers, and 



< E(yi6ia2 - 7) = blEXfl^x-,>eb-^,y^} = ^(l)" 



□ 
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